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Abstract 

We exploit magnetostriction in polydisperse ferrofluids in order to generate 
nonlinear responses, and apply a thermodynamical method to derive the de- 
sired nonlinear magnetic susceptibility. For an ideal gas, this method has 
been demonstrated to be in excellent agreement with a statistical method. In 
the presence of a sinusoidal ac magnetic field, the magnetization of the poly- 
disperse ferrofluid contains higher-order harmonics, which can be extracted 
analytically by using a perturbation approach. We find that the harmonics are 
sensitive to the particle distribution and the degree of field-induced anisotropy 
of the system. In addition, we find that the magnetization is higher in the 
polydisperse system than in the monodisperse one, as also found by a recent 
Monte Carlo simulation. Thus, it seems possible to detect the size distribution 
in a polydisperse ferrofluid by measuring the harmonics of the magnetization 

under the influence of magnetostriction. 
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I. INTRODUCTION 



Ferrofluids (magnetic fluids) are colloidal suspensions containing single domain nano- 
size ferromagnetic particles dispersed in a carrier liquid ffl. These particles are usually 
stabilized against agglomeration by coating them with long-chain molecules (steric sta- 
bilization), or decorating them with charged groups (electrostatic stabilization). Since 
these particles can easily interact via applied magnetic fields, which in turn can affect 
the viscosity and structural properties, ferrofluids possess a wide variety of potential ap- 
plications in many fields ranging from mechanical engineering |]2],|3| to biomedical appli- 
cations @,[|]. Thus, ferrofluids have received much attention in the scientific commu- 

nity @0,|B0ElElElEMlEi0fflli3- 

Polydispersity of ferrofluids emerges naturally since the particles in them always pos- 
sess a log-normal distribution It has been experimentally observed that different 
microstructures can spontaneously form in ferrofluids [[!(]]. This has a strong effect on the 
macroscopic properties of ferrofluids. In this regard, the influence of polydispersity on the 
magnetization of ferrofluids is of both academic and practical interest. Due to the polydis- 
persity of ferrofluids, their structure and magnetization properties may significantly differ 



from those of monodisperse systems |Tl|JT^ , p!3|jr^] . 

The structure of polydisperse ferrofluids has been discussed theoretically on the basis 
of a bidisperse model in which the fluids consist of two fractions of magnetic particles with 
significant size differences • I n a bidisperse ferrofluid, the smaller particles are affected 

by Brownian motion, and are therefore more or less randomly dispersed. The larger magnetic 
dipole moment of the larger particles leads, however, to a strong interparticle force which 
dominates over Brownian motion. Thus the salient structure in these systems is proposed 
to be a chainlike aggregate formed by the larger particles. Some small particles may be 
attached to the ends of these aggregates, but most of them remain in the non-aggregated 



state [[17]]. These features have also been observed in Langevin dynamics simulations [18 



In addition, Wang and Holm |18|] found that the smaller particles hinder the aggregation 
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of larger particles. Since then, that effect has been investigated in detail by Zubarev ||21|| . 
and the influence of polydispersity on the equilibrium properties of ferrofluids was very 
recently investigated by Kristof and Szalai using Monte Carlo simulations JOJ. Kristof and 
Szalai found that magnetization is generally higher in a polydisperse system than in the 
corresponding mono disperse one. 

An inhomogeneous magnetic field H exerts a translational force F on a ferromagnetic 
particle given by 

F = aH VH + m VH, (1) 

where m and a are the permanent magnetic dipole moment and the magnetizability of 
the particle, respectively. Thus, if the permanent moment points in the direction of H, 
the particles will be displaced towards the regions of higher field strength. In a macroscopic 
sample the average moment is in the direction of the field, i. e., the particles favor orientations 
where their permanent magnetic dipole moments are directed along the field. Thus, an 
inhomogeneous field acting on a macroscopic sample causes a concentration gradient with 
high concentrations at high field strengths. If a sample is situated partially in a field and 
held at constant pressure, the particle density in the field region will increase leading to an 
increase of the permeability. This effect is called magnetostriction, or in general, a response 
of the solution to an inhomogeneous magnetic field. Magnetostriction has been extensively 
studied, e. g., in dipolar fluids |p2|| , single crystals of the high-T c cuprate Bi 2 Sr2CaCu 2 8 p3| . 



polycrystalline Fe films ||24|| , and in cylindrical type-II superconductors p5| . Unfortunately, 
to the best of our knowledge, so far there is neither theory nor experiments dealing with 
the important problem of magnetostriction of ferrofluids. The only notable exception is 
ferrogels |26| which are chemically cross-linked polymer networks swollen with a ferrofluid. 

To model ferrofluids we use a thermodynamical method to derive the magnetostriction- 
induced effective third-order nonlinear magnetic susceptibility £. As a sinusoidal ac magnetic 
field is applied, the magnetization will, in general, consist of ac fields at frequencies of the 
higher-order harmonics. We derive the harmonics analytically by using a perturbation ap- 



proach. The aim of the present paper is to exploit magnetostriction in ferrofluids in order to 
generate a nonlinear response and thus a harmonic response. In experiments, measurements 
of the harmonic responses of magnetization have been used to obtain information of the 
anisotropy distribution in a ferromagnetic amorphous alloy |27J . 

This paper is organized as follows. In Sec. II, the nonlinear magnetic susceptibility arising 
from the magnetostriction is derived, and the harmonics of the magnetization are extracted 
analytically. In Sec. Ill, we numerically calculate the harmonics of the magnetization under 
various conditions. We finish with a discussion and conclusions in Sec. IV. 



II. FORMALISM 
A. Nonlinear characteristic arising from magnetostriction 

Let us assume that a ferrofluid is placed in an inhomogeneous magnetic field and kept 
at constant pressure. Then, the density of the particles in the field regions will increase due 
to the interaction between the permanent magnetic moments of the particles and the field 
leading to an increase in the permeability. This effect is called magnetostriction. 

The experimental situation is the following: There is a field-affected volume with volume 
V c , in which the magnetic field and the magnetic induction are denoted by H c and B c , 
respectively. Both of them should satisfy the usual magnetostatic equations PS|] , 



V • B c = 0, (2) 
V x H c = 0. (3) 

Here Eq. (|3|) implies that the magnetic field H c can be expressed as the gradient of a magnetic 
scalar potential $ such that 

H c = -V<£>. (4) 

Under appropriate boundary conditions, the inhomogeneous ferrofluid inside the field- 
affected volume can be represented as a region of volume V c , surrounded by a surface S'. 
Such boundary conditions can be written as 



$ = — H • X on S', (5) 

which, if the ferrofluid within V c were uniform, would give rise to a magnetic field which 
is identical to H (external field) everywhere within V c . In fact, this boundary condition 
guarantees that even in an inhomogeneous ferrofluid the volume average of the magnetic 
field (H c ) within V c still equals that of the external field, (H), namely, 

(Hc) = ^/ Hc(X)d3x = (H) - (6) 

In this case there is no external field outside the field-affected volume and (or, in practice, 
the external field in other areas is weak enough to be neglected). The ferrofluid with volume 
V is situated such that it has regions both inside and outside the field-affected volume at a 
constant pressure p. 

In the presence of an inhomogeneous external magnetic field H, the effective linear 
permeability fi e and the effective third-order nonlinear magnetic susceptibility £ for the 
ferrofluid inside the field-affected volume are defined as 

(B c ) = (/i e + 47r£(H) 2 )(H), (7) 

where (• • •) denotes the volume average. Eq. (|7|) implies that there is a nonlinear relation 
between the magnetization, M, and the magnetic field, (H). This will be explicitly shown 
later in Eq. (|3l|) . Further, it is worth noting that the nonlinear term £ of the magnetic 
susceptibility is actually an effective quantity. It appears due to the fact that the particles 
can receive a translational force which drives them into the field-affected volume in the 
presence of an inhomogeneous field. In this regard, rather than H we should use (H) in 
Eq. ([?D in order to derive the (effective) nonlinear term £. 

Alternatively, based on thermodynamics, the permeability /x# including the incremental 
part due to the magnetostriction can be defined as 

"» = (Hy) Tr = ( w) T , + LI m (Wi) T , (H , (m) T , P id ' (8) 
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where p stands for the density of the part of the ferrofluid inside the field-affected volume, 
and d m i n (or <i max ) denotes the minimum (or maximum) particle diameter. Here the size 
distribution of particles f(d) satisfies the known log- normal law P|29|j 

(9) 

\/2nad 2a z 



where a is the standard deviation of In d and S the median diameter. 

In Eq. (H), ( ) T represents the effective permeability including all nonlinear effects 
at a constant density. As a matter of fact, regarding both Eqs. (|7|) and (§), the total 
effective third-order nonlinear susceptibility generally contains two contributions. The first 
is the magnetostriction-induced one considered in this work, and the other is the normal- 
saturation contribution. At large field intensities, the higher terms of the Langevin function 
should be taken into account. This contribution is negative, and is called normal saturation. 
In contrast, the magnetostriction has a positive contribution. In this work, we assume that 
the field is moderate such that the normal-saturation contribution is weak enough to be 
neglected. It should be noted that the argument of the Langevin function can become large 
for a very small number of large particles in the tail of the size distribution. Further, we 
shall also neglect the normal-saturation contribution resulting from the very small amount 
of the large particles since this contribution might be expected to have a very weak effect 
on the total effective third order susceptibility. To summarize the above, throughout this 
work, only the magnetostriction-induced contribution is considered. 

Using the above assumptions ( ) T = ^ e can ^ e gi ven by the anisotropic Clausius- 



Mossotti equation, namely |30 



= a jl f d - mm ( a{d) + r^n ddj (io) 

where /i2 denotes the permeability of the carrier liquid, mo(d) [or a(d)] the permanent mag- 
netic dipole moment (or magnetizability) of a particle with diameter d, N(d) is the number 
density of particles with diameter d, ks the Boltzmann constant, and T the temperature. 
Regarding Eq. (|10D, more issues should be remarked. It is known that the usual (isotropic) 
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Clausius-Mossotti equation does not include the particle-particle interaction. When Lo 
and Yu studied the field-induced structure transformation in electrorheological solids, they 
succeeded in developing a generalized Clausius-Mossotti equation by introducing a local- 
field factor f3' which reflects the particle-particle interaction between the particles in an 
anisotropic lattice. This generalized Clausius-Mossotti approach [Eq. fllUD] is a self-consistent 
determination of the local field due to a lattice of dipole moments. In other words, Eq. (|ToD 
should be expected to include particle-particle interactions (at least to some extent), and the 
degree of the particle-particle interactions depends on how much gi deviates from 1/3 (note 
that here Ql = f3' /3), where gi represents the demagnetizing factor in the longitudinal field 
case. In particular, the case when 1/ gi = 3 (or g L = 1/3) corresponds to the isotropic case, 
which yields the well-known (isotropic) Clausius-Mossotti equation. It is worth noting that 
there is a sum rule gi + 2 g? = 1 |31j , where gx is the demagnetizing factor in the transverse 



field case. For electrorheological fluids, similar factors were measured in simulations ||32| , |33 



In Eq. (fTtH), the term results from the average contribution of the permanent 

magnetic dipole moment to the average value of the work required to bring a particle with 
diameter d into the field (H) . More precisely, the mean value of the component of the dipole 
moment in the direction of the field is given by 

mo(d)m = 9^<#>, (11) 

with 7 = m °^y H ^ • That is, we set the Langevin function to L(j) = 7/3. Regarding this rela- 
tion, the following issues should be noted. In the present work, we shall adopt a perturbation 
approach |34| [see Sec. 11(C)], which is suitable for a weak nonlinearity. In the perturbation 



approach, it is well established that the effective third-order nonlinear susceptibility can be 



calculated from the linear field ||35|| , while the effective higher-order nonlinearity must de- 
pend on the nonlinear field []34]. In fact, we could have adopted a self-consistent mean-field 
approach ||36|| , but the perturbation approach appears to be more convenient for analytic 



expressions |3TSfl. Thus, to be able to focus on (weak) third-order nonlinearity, it suffices to 



use the Clausius-Mossotti equation [Eq. (|10|)1 by taking into account the linear relation only, 



i.e., L(y) = 7/3. Due to the same reason, in what follows we shall omit the contribution 
from the nonlinear field, too. 

Regarding the incremental permeability due to the magnetostriction, namely the last 
term in Eq. @ 



<H) \"\"/ / T,p 



we obtain 



This equation is valid for the lowest-order perturbation. The differential increase of the 
density inside the field-affected volume dp(d) corresponds to an increase in mass given by 
V c dp(d). This increase is equal to the decrease in mass outside the field-affected volume 
given by — p(d)d(V — V c ) — —p(d)dV, so that dp(d) = — [p(d)/V c ]dV. Thus, we may rewrite 
Eq. (|12p as 

^-t/M»(fe) M f(S)^ 

Next, we can obtain (^gy) T based on the differential of the free energy dF: 

dF = -SdT-pdV + ^(H)d(B c ), (14) 

4:71 

where S denotes the entropy. Introducing the transformed free enthalpy G, 

G = F + pV-^(U}(B c ), (15) 

we obtain the following expression for its differential: 

dG = -SdT + Vdp-— (B c )d(H). (16) 

Atx 



From this equation, we find 

( dV \ V c fd(B c )\ V C (H) fd(i e 



d{U)J Tp 4vr V dp t 

> r,(H) ^ n V ®P J r,(H) 



(17) 



Then, the substitution of this into Eq. (|13D leads to 

12^(H) 2 = f{d)^p{d) [ d -f] dd. (18) 

Jd min Ait '\d P (d)) TAn) \dp ) T{n) 



We now use 



\dp) T{n) \dp(d)) TAH) ( dp ) T{H) {dp(d)) T ' (19) 

where (3 = —-^ (^~) T is the compressibility in the absence of the field. In the last term of 
Eq. flT9|), terms depending on (H) have been neglected since they lead to terms in powers of 
(H) higher than the second term in Eq. ([18]). In the same approximation, after substitution 
of Eq. QTJ) into Eq. ([18]) we obtain 

/ '-Wx /TT\ 2 _ / Fill, \ 2 



Thus, we have 

t^C: m ^(M T id - (2i) 

So far, the effective third-order nonlinear magnetic susceptibility £ has been derived in 
terms of the compressibility, size distribution function, and the differential of effective linear 
permeability with respect to the density. 

B. Comparison with a statistical method 

The increase of the density Ap due to magnetostriction can be calculated from 
(dV/d(U)) Tp . Details follows: 

/-(H) fdp{d)\ 



ld min Jo v c \a(H) 1 



<W r(H) (H)p(d) (dp e 



d(H)dd 



d min Jo 47T \dp J T (H) 

Jd mia Jo 4tt \dp(d)J T w v ' 



To obtain this equation, Eq. ([19]) has been used, which means that in the expression for Ap 
terms in powers of (H) higher than the second have been neglected. Then, the integration 
can be performed, and we obtain 



Ap 



r-f^w^) di (23) 

7d min 8tt \dp{d)) T 



For comparing with a statistical method, let us apply our formalism to a (monodisperse) 
ideal gas For its effective permeability, in view of gi = 1/3, Eq. (|K]) can be rewritten 



as 



p e + 2 3 \ 3k B T 
For the ideal gas, p e ~ 1, and hence we obtain 



^- 1 --N(a + ^-). (24) 



dp e \ 4tt ( m. 



\op J T m m \ 3k B T J 

where m m denotes the mass of a single molecule, m m = M/Na- Here M is the molecular 
weight and Na the Avogadro constant. For this case, using (5 = M/pRT, Eq. ( |2"BD predicts 

A/? = ^^( a+ 3^Tj- (26) 

Here R represents the molar gas constant. 

This result [Eq. ([26])] can also be achieved by using a statistical method. Because of 
Boltzmann's distribution law, , the number of moles per cm 3 of the gas at a point with 
field strength (H), is given by 

n (H> = n (o) exp (_ W // fcs T), (27) 

where is the number of moles per cm 3 of the gas at a point with zero field, and W the 
average value of the work required to bring a molecule into the field (H), 

w -\(° + 4f) {K? - (28) 

Thus, we have 
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= exp 



iL + ^M (H) 2 /k B T 
2 I 3k B T I 



(29) 

Let us neglect the terms in higher than second powers of (H), and we have 



A„=M(n<«> -„«») = ^|^( 



ml 



For an ideal gas, Eq. fl3D| ) yields exactly the same result as Eq. fl2"E|), albeit derived using a 
different approach. This shows the consistency of our arguments. 

C. Magnetization and high-order harmonics 

The orientational magnetization (z axis) has the following general form 

M= ^V^(H) + £(H) 3 . (31) 

47T 

Here, the higher order terms have been omitted. We use an inhomogeneous sinusoidal ac 
field H = (l z /L z )H ac (t)z = (l z /L z )H ac zsm(ut), where < l z < L z with L z being the length 
of the field-affected volume along z axis. Without the loss of generality, we set L z = 1 in 
the following. We now obtain 

M = ^LZJ^H^t) + ^ ac (t) 3 . (32) 

In view of H ac (t) = H ac sm(u>t), the magnetization M can be expressed in terms of the 
odd-order harmonics such that 

M .U_. sin(wt) + M 3uj sin(3u/£), (33) 

where the fundamental and third-order harmonics are given by 

= ^ZJ^l Hac + ^Hi, (34) 
8tt 32 ac v ; 

M*, = - J^L (35) 
In the above derivation, we have used the identity sin 3 (cut) = (3/4) sin(u;t) — (1/4) sm(3ut). 
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III. NUMERICAL RESULTS 



Without any loss of generality, we choose the following parameters for our numerical 
calculations: /i2 = 1 (non- magnetic carrier fluid), density of the bulk material of the particles 
7g/cm 3 , # ac = 20 Oe, (3 = 0.62 x 10" 10 cm-s 2 /g, a(d) = (due to the small size of the 
particles), and T = 298 K. In addition, the volume fraction of the particles is set to be 0.08, 
and the saturation magnetization of the bulk material of the particles is 450 emu. Finally, 
setting d m i n = 1 nm, and d max = 30 nm ensures 



as expected. 

Based on the model parameters, we calculated the dipolar coupling constant |37] A (5) = 
iM?kBT8 3 ' an d found A(9.5 nm) = 1.16, A(10 nm) = 1.35, A(10.5 nm) = 1.56, which ensures the 
assumption that the particle interaction in our system is weak. 

In Fig. 1, we display the size distribution of the particles in the lognormal law for different 
(a) median diameter 6 and (b) standard deviation a. Fig. 2 shows the fundamental [Fig. 2(a)] 
and third-order [Fig. 2(c)] harmonics of the magnetization as a function of the degree of 
anisotropy 1/gi for different median diameters 5. The size distribution of the particles is 
shown in Fig. 1(a). It is found that increasing the degree of anisotropy 1/ gi causes both 
the fundamental and third-order harmonics to increase. Also, a higher median diameter S 
leads to larger harmonics. 

In Fig. 3, the fundamental [Fig. 2(b)] and third-order [Fig. 2(c)] harmonics of the magne- 
tization are plotted as a function oil/ gi for different standard deviations a. The lognormal 
size distribution of the particles is shown in Fig. 1(b). Again, it is shown that the harmonics 
increase with increasing median diameter S. 

Finally, to compare the above polydisperse case with the corresponding monodisperse 
one, we study the monodisperse case in Fig. 4 for three different diameters which have the 
same values as the median diameters used in Fig. 2. In the monodisperse case, it is also 
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evident that increasing the degree of anisotropy 1/gi causes both the fundamental and third- 
order harmonics to increase. In addition, larger diameter leads to larger harmonics. It is 
worthing noting that both the fundamental and third-order harmonics of the magnetization 
are higher in the polydisperse system than in the monodisperse one when comparing Fig. 2 
with Fig. 4. In particular, the third-order harmonics of the polydisperse system can be of 
two orders of magnitude larger than those of the monodisperse system. In other words, the 
magnetization is higher in the polydisperse system than in the monodisperse one, due to 
the fact that for this comparison the volume fraction of the particles is fixed. This is in 
agreement with the findings of Ref. [19| where a Monte Carlo simulation was used to study 
the influence of polydispersity on the equilibrium properties of ferrofluids. 



IV. DISCUSSION AND CONCLUSION 

Here some comments are in order. In the present paper, we have exploited magnetostric- 
tion in ferrofluids in order to generate nonlinear responses. The proposed mechanism should 
work for dc magnetic fields. It will also work for ac fields with frequency v = uj/(2ir) if the 
size of the sample is not greater than c s /u, where c s is the sound velocity. Thus, v can be 
up to kHz or so. Otherwise the required mass density oscillations will not be able to keep 
up with the rapid changes in the magnetic field. 

To obtain the lowest-order (i. e., cubic) nonlinearity, we have assumed that material prop- 
erties such as permeability of the polydisperse system can be calculated as a linear superpo- 
sition of the corresponding values in the monodisperse systems, see Eqs. fllO|) and (p2|). 
For Eq. (^|), the linear superposition should hold since the nonlinear term £ is actually an 
effective quantity which results from all the monodisperse systems. For Eq. fllO|), we used the 
linear superposition again. The reason is that the right-hand side of Eq. actually repre- 
sents the effective contribution from two parts: The induced magnetization (which has been 
assumed to disappear due to the small size of the particles in our numerical calculations), 
and the permanent-moment-related magnetization. In addition, once the inhomogeneous 
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field is applied, the particles with different sizes are able to move into the field-affected 
volume, thus yielding an increasing particle density. In this regard, for Eq. (p^), the linear 
superposition should be used naturally. 

Nonlinear optical materials with large values of (effective) third-order nonlinear dielectric 
susceptibilities are in great need in industrial applications such as nonlinear optical 



switching devices for use in photonics, and real-time coherent optical signal processors, and 
so on. Due to the similarity between magnetics and dielectrics, the present (effective) third- 
order nonlinear magnetic susceptibilities are expected to have some potential applications 
in nonlinear magnetic devices. 

Throughout the paper, only odd-order harmonics are induced to appear. As a matter of 
fact, if one applies an ac magnetic field superimposed to a dc field, the even-order harmonics 
should appear p7|. That is due to the coupling between the two kinds of fields. On the 



other hand, since the second-order harmonics are often of several orders of magnitude larger 
than the corresponding third-order, the second-order harmonics are more attractive for the 
experimental measurements | |27| |. 

We have considered the fundamental and third-order harmonics. In fact, we can consider 
much higher-order harmonics pB|,|5P|, such as fifth-, seventh-, etc. In doing so, we need to 



keep more terms in powers of (H) higher than the second in Eq. (|18|) . Accordingly, more 
terms should be included in Eq. ([H]). However, such higher-order harmonics are often of 
several orders of magnitude smaller than the third-order, and thus less attractive. 

In the numerical calculations, we have omitted the magnetizability of the particles due 
to the fact that the sizes of the particles are very small in ferrofluids. For these particles, 
the permanent magnetic dipole moments play the main role. However, in case of a magne- 
torheological fluid, the magnetizability of the particles should be taken into account since 
the particle sizes range from 2 to 20 fim - about three orders of magnitude larger than in 
ferrofluids. Fortunately, for treating magnetorheological fluids, the present theory holds as 
well. 

In this paper, we have investigated a log-normal distribution [see Eq. @]. Our the- 
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ory could be extended to treat other particle distributions as well. For instance, for a T 



distribution pEOfl , we should replace Eq. (|9]) with 



do \do J T(a + 1) 

where do and a are the parameters of the distribution, and F denotes the gamma function. 

In view of the existing theory which deals with the influence of polydispersity on the 
magnetization of ferrofluids, such as a perturbation theoretical study [|T2] and a cluster 
expansion study it is instructive to compare the present theory with these methods. 

To sum up, we have used a thermodynamical method to derive the nonlinear magnetic 
susceptibility resulting from magnetostriction, which further yields the harmonics of mag- 
netization in response to an applied ac magnetic field. For an ideal gas, this method has 
been shown to be in excellent agreement with a statistical method. It has been shown that 
the harmonics are sensitive to the particle distribution (namely, median diameters and stan- 
dard deviations) and degree of field-induced anisotropy of the system. In addition, we also 
find that the magnetization is higher in the polydisperse system than in the corresponding 
monodisperse one, which is in agreement with previous findings. Thus, it seems possible to 
detect the size distribution in the polydisperse ferrofluids by measuring the harmonics of the 
magnetization of colloidal ferrofluids under the influence of magnetostriction. In detail, the 
size distribution might be achieved by using Eq. (^) and choosing a suitable distribution 
for f(d) to fit experimental data. 
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FIGURES 

FIG. 1. Lognormal distribution of particles for different (a) median diameter 5 and (b) standard 
deviation a. Parameters: (a) a = 0.45 and (b) 5 = 10 nm. 

FIG. 2. Poly disperse case, (a) Fundamental and (b) third-order harmonics of the magnetization 
against the degree of anisotropy 1/gL for different median diameter 5. Parameter: a = 0.45. 

FIG. 3. Polydisperse case, (a) Fundamental and (b) third-order harmonics of the magnetization 
against the degree of anisotropy 1/gL for different standard deviation a. Parameter: 5 = 10 nm. 

FIG. 4. Monodisperse case, (a) Fundamental and (b) third-order harmonics of the magnetiza- 
tion against the degree of anisotropy 1/gL for different diameter d. 
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